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ABSTRACT 

In an attempt to reconcile recent spectral data with predictions of the standard cooling 
flow model, it has been suggested that the metals in the intracluster medium (ICM) 
might be distributed inhomogeneously on small scales. We investigate the possible 
consequences of such a situation within the framework of the cooling flow scenario. 
Using the standard isobaric cooling flow model, we study the ability of such metallic- 
ity variations to preferentially suppress low-temperature line emission in cooling flow 
spectra. We then use simple numerical simulations to investigate the temporal and 
spatial evolution of the ICM when the metals are distributed in such a fashion. Simu- 
lated observations are used to study the constraints real data can place on conditions 
in the ICM. The difficulty of ruling out abundance variations on small spatial scales 
with current observational limits is emphasized. We find that a bimodal distribution 
of metals may give rise to interesting effects in the observed abundance profile, in that 
apparent abundance gradients with central abundance drops and off-centre peaks, sim- 
ilar to those seen recently in some clusters, are produced. Different elements behave 
in different fashion as governed by the temperature dependence of their equivalent 
widths. Our overall conclusion is that, whilst this process alone seems unlikely to be 
able to account for the sharp reduction in low temperature emission lines seen in cur- 
rent spectral data, a contribution at some level is possible and difficult to rule out. The 
possibility of small-scale metallicity variations should be considered when analysing 
high resolution cluster X-ray spectra. 

Key words: cooling flows - galaxies: abundances - galaxies: clusters: general - 
X-rays: galaxies 



1 INTRODUCTION 

Data from the latest generation of X-ray satellites, Chan- 
dra and XMM-Newton, are forcing us to re-examine some of 



FabiarJ 1994). Various authors f 


e.g., Kaastra et al. 


2001; 


Pe- 


terson et al.| 


2001; 


Tamura et al. 


2001) have drawn attention 



to the discrepancy between the predictions of the standard 
cooling flow model and observed spectra of cluster central 
regions. The expected emission lines for several important 
species (e.g. the Fexvn 15 and 17 A lines) do not seem to 
be present at the levels which simple models would expect. 
This is a trend that appears to be repeated in many clus- 
ters th at have traditionally been thought to harbour cooling 
flows (Peterson et al. 2002). Lines such as these are impor- 
tant because they are strong indicators of low (< 1 keV) 
temperature gas (see Section 2.1). In a standard cooling 



flow, in which gas is cooling down to essentially zero, we 
expect a significant flux in such lines. 

Several ideas (e.g., 



Fabian et al. 



2001 



Peterson et al 



2001) have been put forward to explain this discrepancy. 
Here we focus on just one of these, the suggestion that the 
intracluster medium (ICM) metals might be distributed in- 
homogeneously on small- scales. We examine in detail the 
consequences of such a scenario, investigating the effects of 
an ICM metallicity which varies on small, unresolved scales 
(sub-kpc, say). This idea is no more than a minor extension 
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of the multiphase cooling flow model (Nulsen 1986), which 
has always relied upon the coexistence of phases with dif- 
ferent densities and temperatures at the same radius in the 
ICM. We merely allow the chemical composition of these 
phases to vary as well. This idea is attractive in the context 
of cooling flows since it leaves the spectra of high (> 1 keV) 
temperature gas essentially unchanged, whilst reducing the 
line emission from low temperature gas. It (potentially) en- 
ables us to reconcile the data and models whilst making the 
smallest conceptual changes to the standard model. 
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Indeed, given that a suppression of thermal conduction 
(electron motion) implies an even more severe reduction 
in the freedom of heavy ions to diffuse through the ICM, 
we might expect small-scale metallicity variations to be a 
natural consequence of the multiphase cooling flow model, 
in which conduction is by necessity heavily reduced. Con- 
versely, if ICM thermal conduction is relatively e fficient (as 
has bee n suggested recently by sever al authors, e.g., Narayan 



Med 



edev 



[2001 ; Voigt et al. 2002) then this does not nec- 



essarily imply that the same is true of ion motion. 

In Section 2 we examine the effects of an inhomogeneous 
metallicity within the fr amework of the stan dard isobaric 
cooling flow model (e.g., Johnstone et al, 1992), quantifying 



the reduction in equivalent width that can be expected for 
several important ICM species as the degree of inhomogene- 
ity is varied. In Section 3 we outline an improved numerical 
model that allows us to carry out simulations with spatial 
and temporal evolution, and to produce simulated obser- 
vations of Chandra spectral data. Section 4 describes some 
of the results of this model for a simple case of metallicity 
inhomogeneities. We find that small-scale abundance varia- 
tions may give rise to the appearance of abundance gradi- 
ents where none in fact exist, with the results for individual 
species being controlled by the temperature dependence of 
their equivalent widths. We stress the difficulty of resolving 
such variations without a combination of good spatial and 
spectral resolution. In the final Section we make some dis- 
cussion concerning the likelihood of such small-scale metal- 
licity variations in the ICM. Given that direct detection of 
such structures is currently unfeasible, we put forward some 
suggestive circumstantial evidence. 

Throughout this paper we assume: Ho = 
50kms _1 Mpc" 1 ; Qm = 1.0, = 0.0. We make use 
of XSPEC version 11. 
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Figure 1. Behaviour of the Fexvn 15.0 A line in a coronal 
equilibrium single temperature plasma. The curve labels refer to 
metallicity in solar units. The top panel shows the equivalent 
width in eV, the bottom panel shows the fractional power (in the 
X-ray waveband 0.1-10.0 keV) contained in this line. 



2 COOLING FLOW SPECTRA 
2.1 ICM thermometers 

Fig. [l^, shows the temperature dependence of the equiva- 
lent width of the 15 A Fexvn line, for various metallicities, 
for a single temperat ure plasma, as ca l culated using the 



meka l plasma code (Mewe et al 



19!).-) 



Kaastra & Mewe 



1992 ) . This code computes the bremsstrahlung continuum 
and line radiation from a plasma (of up to 1 5 elements) 
in the coron al equilibrium approximation (e.g., Brickhousc 

10 a K 



et al 



The 



1995) over the temperature range 10 4 
equivalent width curves are sharply peaked around 0.3 keV 
(~ 10 6 K), so these emission lines are sensitive indicators of 
gas at such temperatures. Furthermore, they are extremely 
strong - Fig. |l|b shows the fraction of the power (in the 
X-ray waveband 0.1-10.0 keV) being emitted at any tem- 
perature in this iron line. For a solar metallicity plasma, it 
peaks around the 10 per cent level. The results are similar 
for the 17 A line. Consequently, these lines are highly im- 
portant coolants below 1 keV, and therefore in the standard 
cooling flow model (in which gas cools from a temperature 
of several keV to essentially zero) , we expect a strong flux in 



ters are much weaker than would be expected ( 


3.g., [Kaastra 


et al. 2001 




Peterson et al. 


2001 




Tamura et al. 


2001). 



2.2 Bimodal metallicities 

2.2.1 Fixed temperature plasn 
Several ideas (e.g. 



Fabian et al 



2001 



Peterson et al 



2001) 



have been put forward to explain this missing low tempera- 
ture flux. Here we focus on just one of these, the suggestion 
that the ICM metals might be inhomogeneously distributed 
on small-scales. The essence of this idea is as follows. At any 
fixed single temperature, a two-component plasma consist- 
ing of metal-rich and metal-poor 'phases' (note that the us- 
age of the term phase here is not necessarily strictly the same 
as in traditional multiphase cooling flow models) is spectro- 
scopically indistinguishable from a homogeneous plasma of 
some mean metallicity Z — /hiZhi + (1 — fhi)Z\ , where the 
metal-rich phase has a metallicity Zhi and accounts for a 
mass fraction /hi, and the metal-poor phase has a metallic- 
ity Z\ a . For simplicity (and without loss of generality), we 
take Z\ = 0. /hi = 1 corresponds to a uniform plasma, and 
decreasing values of /hi to increasingly segregated plasmas. 
Throughout this paper, we leave the He abundance fixed at 
the solar value, allowing the heavy elements to va ry as spec- 
ified by Z. We use the solar abundance ratios of 
Grevesse (1989), hereinafter Zq. 



Anders & 



The indistinguishability of the spectra in these two cases 
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is a simple consequence of the fact that the strength of the 
continuum radiation is independent of Z, whilst that of the 
emission lines is directly proportional to Z. The reduction 
in the mass fraction of the enriched phase is therefore offset 
by the increased line strength. Alternatively, in the coronal 
limit, the local concentration of the heavy elements makes no 
difference to their radiation: a given number of heavy ions 
will radiate in the same way whether uniformly dispersed 
throughout the emitting volume or concentrated in a specific 
region. As an aside, we note that this simplistic treatment 
breaks down at very high metallicities (Z > 70 Zq) when 
the free electrons contributed by the metal ions become sig- 
nificant. This would not invalidate the idea that a mixture 
of a metal-rich and a metal-poor plasma can look identical 
to a uniform metallicity plasma, it would merely change the 
dependence of the normalization of each component on Z 
from a simple linear relationship to something more compli- 
cated. Such ultra-high metallicities are of no relevance here 
though. 



2.2.2 Cooling plasmas 

The situation becomes more complex when one allows cool- 
ing to occur. We briefly review th e derivation of the sta ndard 
isobaric cooling flow model (e.g., Johnstone et al. 1992). The 
energy released per unit mass of gas on cooling by d T is 



dq 



d u + pd I — 



3 k B T 



+ 



= d 
nk B T 



u + - 



(constant p) 



(ideal gas). 



^ 2 fimn p 
For a constant mass-flow rate M , the power release is 



thus 
dL = l 



fcp 



-MdT. (1) 

2 /ifflH 

We also have, from the definition of the cooling function 



A. 



dL = n c nnA(T,Z)dV, 
dL u = n e rmKu{T,Z)dV. 



(2) 



Hence we obtain the spectral power for a steady-state 
flow cooling from T max to T m i n 



-— 2_M 

2^m H JT 



A„(T, Z) 
A(T, Z) 



dr. 



(3) 



Thus, we have the simple result that the emission mea- 
sure of each temperature component in the flow is inversely 
proportional to the cooling function A(T) at that tempera- 
ture. This model is implemented by the xspec package as 
the mkcf low and vmcf low models^J We have constructed an 
independent implementation of Eqn. [], calculating cooling 



1 We draw attention to the fact that in the XSPEC mkcf low imple- 
mentation, altering the plasma bulk metallicity varies the helium 
abundance (for XSPEC versions prior to ll.l.Oab). This can have a 
non-negligible effect on the calculation of equivalent widths when 
high metallicities are involved. As a workaround, the vmcflow 
model may be used to specify abundances on an element-by- 
element basis. 




Mass fraction of enriched phase (f hi ) 

Figure 2. Fractional change in the equivalent width of some 
important emission lines from the cooling flow model of Eqn. 
for two-component plasmas, as the mass fraction of the enriched 
phase is varied. The dashed line is the reference line y = x. Param- 
eters of the cooling flow model: Z = 0.3 Zq , Z\ Q = 0, Z^i = f^ 1 Z. 
The integration over temperature was performed at 100 points, 
logarithmically spaced between 0.0808-8.620 keV. 



functions by integrating mekal spectra against energy over 
the 5eV-200keV range. 

The heating solutions to the cooling flow spectral prob- 
lem (i.e. the absence of spectral lines characteristic of cool, 
< 1 keV, gas) suppress the low-temperature line emission 
by raising T m i n to ~ 1 keV. The bimodal metallicity hy- 
pothesis suppresses the equivalent width (EW) of the low- 
temperature lines by a different means. In Fig. ^| we illus- 
trate how the EW of various important spectral features 
changes as a function of fhi (or Zhi), for the case Z — 0.3 Zq. 
The dashed line shows what might be the naive expectation, 
namely that EW(/m) = /hi EW(1). In no case is this line 
followed - the reduction in the EW is always less than this. 
Thus, for example, at /m = 0.1, the EW of the Fe Kq line 
is reduced by a factor of about 0.5, that of the 19 A O line 
by about 0.25, and that of the 15 A Fe line by around 0.15. 
We can interpret these results as follows. 

Recall the definition of equivalent width 



EW 



El 



B 



+C 



Ic 



la 



dE, 



(4) 



where Jl+c is the total (line + continuum) intensity at a 
given energy, and 7c is the continuum intensity at the same 
point. 

For fixed, narrow line profiles, we can neglect the in- 
tegral over energy and simply compare intensity ratios at 
the line energy. Making use of the cooling flow equation 
Eqn. [], we may write the equivalent width, EW(H), for a 
given spectral line (with a narrow energy profile centred on 
some energy E) from a homogeneous metallicity cooling flow 
as 



EW(H) = 



A L (T,g,Z) 
A(T, Z) 



dT 



Ac(T,E,Z) 
A(T, Z) 



(5) 



dT 
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In this expression, AL?i c nH would be the power per unit 
volume per unit energy radiated by the line component at 
a given energy, with Ac having an analogous meaning for 
the continuum component. Integrating the sum of these two 
functions over energy leads to the cooling function A(T, Z). 

For a two-component ('bimodal') cooling flow, in which 
one component has a metallicity Zu and accounts for a frac- 
tion /hi of the total mass flow rate, and the other component 
has zero metals, we have the following expression for the 
equivalent width, EW(B), of a narrow spectral line centred 
at energy E: 



EW(B)= / hi 



A L (r,£,z hi ) 

A(T, Zu) 



dT I 



(6) 



/hi 



A C (T,E,Z M ) 
A(T, Zu) 



+ (1 - /w) 



A c (T,£,0) 



A(T,0) 



dT 



Our desire is to suppress the equivalent width of the 
low-temperature lines in the bimodal case - that is, the ratio 
of Eqn. g to Eqn. ^ should be < 1 in such cases. There exists 
no simple algebraic solution to this form, owing to the non- 
analytic behaviour of the A functions. 

Restricting our attention to an emission line that ex- 
ists only over a narrow temperature range (for example, the 
Fexvn 15 A line, as illustrated in Fig. Jp) , however, we may 
make the simplification J d T — > £^ ST. In addition, 
we have the result that AL,n e riH oc n e ni on Al oc m on tx Z 



(7) 



(at fixed tih). Consequently 

A L (f,E,Zhi) = Zin = J_ 

A h {T,E,Z) Z /hi' 

The ratio of the equivalent widths may therefore be some- 
what simplified to 

EW(B) _ A(f,Z) f Tmax A C (T,E,Z ) 
EW(H) A(f,Zhi)7 Tmin A(T,Z) 



-dT / 



(8) 



/hi 



A c (T,E,Z hi ) 
A(T, Zu) 



/hi) 



A c (T,E,0) 



A(T,0) 



dT 



The term involving the ratio of integrals is the ratio of 
the continuum in the homogeneous case to that in the bi- 
modal case. By inspection of cooling flow spectra, this ra- 
tio is very close to one. Most of the continuum flux is con- 
tributed by gas at temperatures > 1 keV, where the metals 
have relatively little influence on the cooling function. Bi- 
modal plasmas tend to have fractionally stronger continua at 
low (< 1 keV) energies, where low-temperature gas makes a 
slight contribution. The reduced cooling function of the low- 
metallicity phase at low temperatures enhances its emission 
measure compared to the mean metallicity case. 

In the idealized case of a narrow emission line, the ratio 
of the equivalent widths for the two cases therefore reduces 
to the ratio of the cooling functions A(T, Z)/A(T, Zu) - For a 
high-temperature (> 1 keV) line, where cooling is essentially 
independent of metallicity, the ratio of the equivalent widths 
reduces to one. For a low-temperature line, then if the metals 
were to completely dominate the cooling, we expect the ratio 
to reduce to Z /Zu = fu- 

Qualitatively, on increasing the metallicity by a fac- 
tor /h^ 1 , the emission line strength (at all temperatures) 
increases proportionately. Reducing the mass fraction com- 
pensates, so that spectra at any fixed temperature are un- 




Temperature / K 

Figure 3. Dependence of the cooling function A(T) on metal- 
licity. The labelled curves show the cooling function for a given 
metallicity (in solar units) relative to that for 0.3 Zq. The num- 
bers in parentheses are the metallicities as multiples of 0.3 Zq. 
At low temperatures, the curves asymptote out to these values. 
Above 10 8 K the cooling is dominated by the H, He continuum; 
as the temperature decreases to ~ 10 6 K Fe becomes increas- 
ingly dominant; and towards the end of the range O is the ma- 
jor coolant. The dashed curve shows the equivalent width of the 
Fexvn 15 A line in arbitrary units. 



changed. When steady-state cooling is allowed to take place, 
however, the faster cooling of the metal-rich gas at low tem- 
peratures reduces its emission measure and suppresses its 
spectral contribution. 

The cooling function ratio is plotted for some represen- 
tative metallicities in Fig. |^. Note that metals do not fully 
dominate the cooling function until T < 10 6 K. In the tem- 
perature range where the Fexvn 15 A line is strong (over- 
laid), the ratio of the cooling function for Z = 3.0 Zq to that 
for Z = 0.3 Zq is in the range 6-7. This is consistent with 
the result from Fig. ^| for this line with fu = 0.1, where a 
fractional suppression of the EW of order 0.15 is found. This 
is to be contrasted with the simple expectation of 0.1. The 
actual suppression of this line is not so strong. 

In practice for lines existing at intermediate temper- 
atures, or over a non- negligible range of temperatures, we 
expect behaviour between these limits. The iron K line, for 
example, deriving as it does fromgas over a relatively wide, 
high temperature range (see Fig. p|), has its equivalent width 
suppressed by about a factor of 0.5, for fu ~ 0.1. If this were 
the case, then analyses not taking this into account would 
have underestimated the mean ICM metallicity. As another 
consequence of this, the relative suppression of the Fexvn 
line (say) compared to that of the Fe K line, is only a fac- 
tor of 0.3. In other words, given a spectrum of unknown 
metallicity, if one were to fit an abundance from the iron 
K line, then attempt to reproduce the observed Fe xvn line 
with this abundance, the observed line would be about a 
factor of three smaller than predicted. This is of course ap- 
preciably less than the simple 'factor of ten' reduction that 
might have been expected. Consequently, it seems unlikely 
that small-scale metallicity variations alone can be responsi- 
ble for the very strong suppression of the equivalent widths 
of low-temperature emission lines seen in cooling flow spec- 
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tra, unless the range of the variations were to be pushed to 
extreme levels. 

/hi = 0.1 corresponds to /hi ~ 3 Zq. This may seem like 
an excessively high metallicity, but consider that supernovae 
can be viewed as sources of essentially infinite metallicity. 
The Chandra observation (Iwasawa et al. |2001 ) of 4C55.15 
reveals a metallicity ~ 2 Zq within the central 50 kpc, as 
well as a metal-rich 'plume' extending over ~ 25 kpc with 
Z ~ 8 Zq . Whilst this is clearly an extreme case, if such 
high-metallicity features can be produced on these scales, 
there is no reason they cannot occur on the much smaller 
scales we consider. 



3 AN IMPROVED COOLING FLOW MODEL 

The standard isobaric cooling-flow model represented by 
Eqn. |3] is a useful tool, but it suffers from a number of lim- 
itations. It assumes a constant pressure, and is inherently 
single-phase. It assumes a steady-state with no time depen- 
dence. Furthermore, there is no spatial information, and no 
attempt to treat the gravitational potential of the cluster. 
It is a good first approximation to the global properties of 
a cooling flow, but in the new era of high spatial resolution 
X-ray spectroscopy and imaging, being able to study the 
temporal and spatial evolution of the ICM becomes more 
and more desirable. An improved model may be o b taine d 
by following the numerical prescription of Thomas ( 1988 ). 
The underlying theoretical basis is that of fNfulsen (1986). 
We assume spherical symmetry, so that the system can be 
described by a one-dimensional model. 



3.1 The dark matter 

We simulate a cluster simply as a two-component system, 
comprised of hot gas in the potential well of a dark matter 
(DM) halo. We represent the latte r by a standard Nav arro, 



Frenk and White (NFW) profile (Navarro et al. 1997). In 



order to specify the parameters of the p rofile, we mak e use o f 
the calibrated virial scaling relations of Evrard et al. (1996): 



(9) 



M(r 200 ) = 2.9 x lO 15 ^, 1 (jq^j) 2 Ms, 

r 200 (T) = 3 .7fc 6 - 1 ( I ^) 5 Mpc; 

where r2oo is the radius at which the density contrast (i.e. 
density in units of the critical density) is equal to 200, and 
serves to defin e a virial radius for the cluster (e.g., Cole 
& Laccy 1996j ). We complete the specification of the DM 
profile by t aking a typical cluster value of 5 (e.g., IBullock 



et al 



2001) for the concentration parameter (i.e. the virial 



radius in units of the NFW scale radius). 



3.2 The gas 

The gas is modelled as a single-phase medium. The tradi- 
tional assumption for cooling flows has been that the un- 
stable nature of the cooli ng process w ould lead to a highly 



multiphase medium (e.g., Allen 2000), i.e. phases of differ- 



Recent observatio ns cast doubt on this assum ption, though 
- see for example Molendi & Pizzolatc (2001). 

For initial conditions, we assume a continuous smooth 
pressure profile of the isothermal form, normalized in terms 
of an (initially arbitrary) outer boundary pressure at the 
virial radius. The gas is assigned a density profile which 
is proportional to the pressure profile. The scaling is such 
that the temperature of the gas at the outer boundary is 
equal to the virial temperature of the cluster. Note that 
when we allow the metallicity of the gas to vary spatially, 
the molecular weight n is not constant, and we calculate 
the outer boundary temperature using a 'mean' /i. Also, in 
such cases the gas will not be isothermal, since for a smooth 
pressure profile to exist, regions of high fi must be at higher 
temperatures. 

The total mass of gas is calculated via numerical inte- 
gration of the density profile. We adjust the outer boundary 
pressure so that the system gas mass corresponds to some 
appropriate fraction / gas of the total (gas plus DM) mass. 
A temperature profile T(r) may then be assigned using the 
(ideal gas) equation of state. 



3.3 Evolution of the system 

The system is allowed to evolve under the influence of grav- 
ity (neglecting the gas self-gravity) and radiative energy loss. 
We simulate the cluster over the radial range < r < r^oo, 
i.e. out to the virial radius. For boundary conditions, we im- 
pose a stationary inner boundary and a constant pressure 
outer boundary. The latter is acceptable because unphysi- 
cally long time-scales would be required for the effects of 
cooling to propagate out as far as the virial radius. The 
governing differential equation s are rep r esent ed by the La- 
grangian difference scheme of Thomas' ( 198S ) . Put simply, 



the radial range is divided into a large number of zones, 
which are then evolved forward in discrete time-steps. We 
use an adaptive step-size that is limited by stability consid- 
erations, but also by the cooling time of the gas, which in 
the standard isobaric assumption is given by 



fc 



( fa 



7-1 V ^ Xh 



TAT 

KTy 



(10) 



ent densities and temperatures coexist at the same radius. 



where: 7 is the adiabatic index of the plasma; Ah = numn/p 
is the hydrogen mass fraction; e r = n c /nu is the elec- 
tron density relative to the hydrogen density. With the 
cooling function A(T) ~ T2 (at high temperatures) for 
bremsstrahlung, then at a fixed pressure t coo i ~ T 2 ; whereas 

for virially-scaled clusters with p oc T at a fixed overdensity, 
'1 

tcool ~ T 2 • 

Since we are concerned only with modelling the gas in 
the X-ray regime, we do not follow the cooling below 10 K. 
At this temperature, the gas is no longer radiating apprecia- 
bly in the X-ray waveband (see Fig. 0), and may therefore 
(for our purposes) be ignored from further consideration. 
In addition, the cooling time (see Fig. ^) is becoming pro- 
hibitively short for continued computation. 

If at the end of an iteration step any zone has cooled 
below 10 5 K, it is removed from subsequent analysis, and 
the remaining zones are adiabatically expanded to fill the 
vacated volume. Zones closest to the loss-site undergo the 
greatest degree of expansion. Conceptually, this is intended 
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Figure 4. Cooling functions of plasmas of varying metallicities. 
He is at solar abundances in all cases. Solid — H, He only; dashed - 
metals at 0.3 solar abundance; dotted - metals at 3.0 solar abun- 
dance. The upper curve in each case is for all thermodynamically 
relevant radiation (5 eV-200 keV); the lower curve is for the X-ray 
wave-band 0.1-10 keV. 



as a (very) crude representation of the rapid collapse of cold 
gas into condensed objects of negligible volume (neglecting 
any kind of energy injection at this point from star formation 
say). 

3.4 The plasma radiation 

We model the radiation coming from the plasma with the 
mekal spectral code. Thus, the abundance of the plasma be- 
comes a relevant parameter, and each radial zone may have 
a different set of abundances. At the start of a simulation, 
spectra for each abundance set are integrated in order to ob- 
tain evaluations of the cooling function A(T), and the cool- 
ing time integral J T at discrete temperatures over 
the range of interest. The relative electron density is also 
calculated as a function of temperature. During the course 
of a simulation we interpolate amongst the tabulated val- 
ues of these functions to find the value appropriate for any 
particular plasma temperature. 

Besides these essential quantities, it is clearly of great 
interest to use mekal to produce synthetic spectra for the 
plasma. These calculations need not be performed at every 
iteration, but only when the physical state of the system is 
to be examined. Note that the energy range of interest is 
different in the two cases. In the case of A(T), we are con- 
cerned with all the radiation that acts as a significant energy 
loss for the plasma over the temperature range of interest. 
In practice we use the energy range 5eV-200keV. In the 
case of spectra, only the much more limited wave-band of, 
say, 0.1-10 keV, is of relevance for comparison with the X- 
ray observations of Chandra and XMM-Newton. These two 
functions are displayed in Fig. ^. Rather than being inter- 
polated, spectra are calculated directly for each zone. 

To facilitate comparison with observations we first in- 
tegrate along a line-of-sight to obtain the 'spectral surface 
brightness' at projected distance b from the centre 

e„(r)2r dr 



Figure 5. Isobaric cooling times for the plasmas of Fig. ^, calcu- 
lated by integrating the total radiation curves for each abundance 
set. The pressure is 10 6 K cm -3 , with the times inversely propor- 
tional to pressure. 



trum due to the (distance-weighted) contributions of all the 
gas along the lines-of-sight between projected radii b 1 and 
b 2 



d b2n b r^± (12) 
lb! Jb Vr 2 -b 2 

47t[F(6i) - F(fta)], where 
dr e v (r)ryj r 2 — b 2 . 



P v (bi O O2 



F(b) 



3.5 The simulated observations 

The integrated spectrum corresponding to each annulus is 
then redshifted in terms of both magnitude (as per the 
appropriate luminosity distance) and energy. The simula- 
tions reported here are for a redshift of 0.017, correspond- 
ing to a luminosity distance of 100 Mpc. These data are then 



conv erted to xspec table-model format fits files (Axnaud 



199E). For added verisimilitude we include the effects of ab- 



(11) 



Vr 2 - b 2 

We then integrate again over an annulus to obtain the spec- 



sorption by multiplying our model spectra with the xspec 
phabs model, using a column density Ah = 10 21 cm -2 . 
Since we are now interested in spatial effects, for response 
matrices we have used the Chandra ACIS-S data, specifically 
the AO-2 proposal planning rmf and ARF filesn for the de- 
tector aim-point position. 

The simulated integration time was typically 25 ks. 
Noise was added via counting statistics, but no background 
file was included. The resulting fake data were grouped to 
a minimum of 20 counts per channel to ensure applicability 
of the x 2 statistic. 



4 A CHEMICALLY INHOMOGENEOUS ICM 

4.1 A bimodal metallicity distribution 

Our aim is to investigate the effects of a 'clumpy' metal 
distribution in the ICM. In our model, each radial zone cor- 
responding to the difference equations may have a distinct 

2 http : //asc .harvard. edu/caldb/Aef f / 
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Figure 6. Time evolution of the observed abundance profile for a bimodal metallicity. Vertical error bars are la. See text for model 
properties. Filled circles - 0.0 Gyr; open squares - 1.0 Gyr; open circles - 1.2 Gyr; open diamonds — 1.6 Gyr. Unless otherwise stated, all 
results are for single temperature models fitted to the 0.3-7.0 keV spectral range. The simulated observation time was 25 ks for the upper 
two panels, and 50 ks for the lower two. Panel (a) : solid lines - single temperature mekal fits; dotted lines — results of two temperature 
fits, shown where the reduced x 2 f° r * ne single temperature fits exceeded 2.0. Panel (b): mekal fits to the 3.0-7.0 keV (i.e. iron K) 
spectral region. Panels (c) and (d): vmekal fits for the silicon and oxygen abundance, respectively. 



set of element abundances. For this work, we have chosen 
one of the simplest possible inhomogeneities: 9 out of every 
10 zones are taken to be pure H, He in the solar abundance 
ratio; whilst every tenth zone has a three times solar abun- 
dance of the 13 heavy elements included by mekal (with 
He remaining fixed at the solar value). A spatially varying 
distribution of this simple form appears identical to a ho- 
mogeneous plasma of abundance a little less than one third 
solar. 

Here we present the results for one particular clus- 
ter, with parameters as follows: virial temperature 8.6 keV 
(10 s K), which from Eqn. |^ corresponds to a virial radius of 
3.4 Mpc and a halo mass of 2.3 x 10 15 Mq. The radial range 
was modelled with 5000 zones, corresponding to a resolution 
of 700 pc, or around 1.5 arcsec for the redshift in qu estion. 
The gas fraction was 0.17 (e.g., |Ettori fc Fabian|[l999[). With 



a concentration of 5, the NFW scale radius is 700 kpc, and 
the associated (singular isothermal sphere) velocity disper- 
sion is 1300kms _1 . 

Fig. |i| illustrates the results of various xspec fits to 
the simulated observations produced for this system. We 
display only the central regions where cooling has progressed 



significantly in the elapsed time. In each case the specified 
XSPEC model was multiplied with a phabs component. The 
redshift z and column density Nh were fixed at the correct 
values. 

Panel (a) illustrates the results for single temperature 
mekal fits to the spectral range 0.3-7 keV, in which the 
plasma temperature, abundance and normalization were al- 
lowed to vary as free parameters. 

Firstly, note that the initial fitted abundance profile is 
fiat (within the given error bars), with no high metallicity 
spikes. This confirms the basic premise that a varying abun- 
dance plasma can resemble a homogeneous gas. Fig. M shows 
the evolution of the true (not fitted) temperature profile of 
the gas with time. Note that in the initial conditions the 
high metallicity zones are at higher temperatures. This is 
a consequence of their higher molecular weight (0.63 as op- 
posed to 0.61 in this specific case). Recall that we require a 
smooth pressure profile, and that this is a function of T/jj,. 
These initially high temperatures are of no consequence for 
the subsequent evolution of the system. 

As the system evolves, it is obviously apparent that the 
metal-rich regions cool more swiftly than the metal-poor re- 
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Figure 7. The physical temperature profiles for the model of 
Fig. ^| In order of decreasing average temperature, the curves 
are: 0.0, 1.0, 1.2 (dotted), 1.6 Gyr. 



gions, leading to the rapid inversion of the 'spikes' in the 
temperature profile shown in Fig. Q This is an inevitable 
consequence of the shorter cooling time of the high abun- 
dance gas (Fig. 5), which in turn is due to its enhanced 
radiation (Fig. U). The cooling is naturally most extreme 
in the central regions where the gas density is highest, and 
the two-body bremsstrahlung radiation is most intense. The 
disappearance of the innermost spikes in the temperature 
profile at late times is due to the removal of the metal-rich 
gas from the simulation as it cools below 10 K. 

This differential cooling has marked effects upon the 
measured abundance profile (Fig. ^a). As time passes, there 
is a gradual increase in the derived abundances, particularly 
in the central regions (for example, after 1.0 Gyr the mea- 
sured value for the central abundance has risen to 0.5 Zq). 
With further evolution, the central abundance begins to de- 
crease again, leading to the situation where there is a peak 
in the abundance profile at an off-centre position. This peak 
declines in magnitude and moves outward with subsequent 
evolution. 

This behaviour may be explained as follows. In our 
simulations, the metal-rich gas (contributing the emission 
lines) cools relatively quickly and enters the regime where 
line emission is more important, whereas the metal-poor gas 
(contributing the bulk of the continuum) cools much more 
slowly. The combined consequence of these two processes is 
an increase in the strength of the lines relative to the con- 
tinuum, that is, an increase in the equivalent width of the 
lines. It is from the equivalent width that the fitting proce- 
dure obtains the plasma abundance. One could view this as 
a form of emission-weighting - as the metal-rich regions cool 
the intensity of their emission increases, so they increasingly 
dominate the spectral fits. The decrease in the central abun- 
dance at late times occurs as metal-rich gas there cools out 
of the X-ray regime and is lost from sight. 

The maximum amplitude of the abundance peak is ob- 
viously related to the underlying spread in metallicity, but 
it is not obvious how one might use the former in practice 
to reconstruct the latter. From the theoretical point of view 
there is a lower limit to the extent of the variations that 
must exist if they are to be of any use in resolving the cool- 
ing flow problem (i.e. the lack of low-energy spectral lines) . 
The upper limit to the spread depends on what one consid- 



ers reasonable for the injection process to produce, assuming 
there is no active segregation of metals. Cooling time con- 
siderations also become important with extreme abundance 
values, if one does not wish the entire central region to be- 
come devoid of metals in the high temperature phases. 

Fig. ^a also displays the results of two-temperature 
mekal fits to the data (dotted lines). These are shown only in 
those cases where the reduced \ 2 f° r the one-temperature 
fits was greater than 2.0. It was necessary to restrict the 
abundances of the two mekal components to be equal, there 
not being enough information in the spectra to obtain mean- 
ingful constraints on two individually varying abundances. 
Even taking this step, the values obtained for the two tem- 
peratures were more often than not very poor limits. The 
resulting abundances, however, agree very well with the one- 
temperature results. 

To check that these effects are a genuine consequence 
of the non-uniform metal distribution, we have carried out 
identical simulations where all the gas has a uniform abun- 
dance of 0.3 Zq. In this case, all four abundance profiles 
analogous to those plotted in Fig. [] remain resolutely flat 
within the errors at all times as the system evolves. 

The results for lower temperature clusters (3 and 5 keV) 
are qualitatively similar. Obviously, with lower temperature 
systems there is a greater need for two-temperature models 
to fit the correspondingly richer spectra. Also, while for the 
8keV cluster the maximum height of the abundance peak 
(0.9 Zq) is some three times the average abundance, for the 
3 keV system the peak has an amplitude roughly double the 
average abundance. There is less opportunity for differential 
cooling if the bulk of the gas begins life at lower tempera- 
tures. The evolution of lower temperature systems also takes 
place more rapidly (recall that from Eqn. [u] that t coo \ ~ T2 
for virially-scaled clusters). The general trends of behaviour 
are, however, the same. The small-scale metallicity variation 
scenario would therefore predict that the off-centre abun- 
dance peaks would be proportionately stronger in higher 
temperature systems. 



4-1.1 Matters of resolution 

The size of the spectral annuli is controlled by two factors. 
Firstly, it is of course necessary to ensure that they are large 
enough to collect sufficient photons in a reasonable integra- 
tion time to produce a meaningful signal-to-noise level. Sec- 
ondly, there is also an issue of numerical resolution. The size 
of the annuli is such that each encompasses ten radial zones 
from inner to outer edge (although of course owing to projec- 
tion effects each annulus receives emission from every zone 
whose radius is greater than that of the inner edge of the 
annulus). This is the minimum number necessary to avoid 
artificial oscillations in the fitted abundance profile, given 
that the fraction of metal-rich zones is 10 per cent. With a 
fraction /hi of zones being metal rich, spectral annuli must 
encompass at least /hi 1 radial zones. 

Alternatively, if the spectral annuli were any smaller 
then we would be able to resolve the discrete nature of the 
metal distribution. Given that each annulus spans a pro- 
jected radius of 7 kpc, whilst the size of the metal rich zones 
is 0.7 kpc, it might seem that we are claiming it would be 
possible to resolve the inhomogeneities with a resolution 
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length greater than their size. This is not the case, how- 
ever, because the spatial extent of the metal-poor zones is 
nine times that of the metal-rich zones. When the resolution 
length is less than or equal to the greater of those extents, 
the distribution can in principle be resolved. 

As is standard practice with numerical simulations, we 
have doubled the numerical resolution to check that the re- 
sults are unaffected, and indeed they are not. This is im- 
portant in this case for another reason. As well as confirm- 
ing the numerical result, it illustrates just how difficult it 
will be to probe such metallicity distributions in real-world 
observations. Looked at another way, then going from the 
high resolution simulation to the default resolution simula- 
tion doubles the physical size scale of the metallicity varia- 
tions, so that the largest relevant length scale (that of the 
metal-poor regions in our model) is just less than the reso- 
lution l ength of the spectral observations. Yet it does not do 



eral clu sters with ASCA and ROSAT, for exampl e: Cen- 
taurus (Fukazawa et al. f994; Alle n fc Fabian | l994 ); Virgo 



anythir g towards hinting that there are in fact extreme en- 
richment variations on just slightly smaller scales. In other 
words, detection of such metal variations will be an all-or- 
nothing affair - without adequate spatial resolution, there 
is no hope. 

Recall as well that this is for a highly simplistic, reg- 
ular, one-dimensional distribution. Given an irregular, two- 
or three- dimensional distribution pattern, with metal-rich 
'clouds' drifting in and out of various lines of sight, then the 
problem becomes even more difficult. 



If ve reduce the size of the spectral annuli below the 



level we have used, then we do not immediately begin to re- 
solve the correct (i.e 0.0-3.0 Zq) metal distribution. Instead, 
one begins to detect only slight variations around the mean 
level. These will not be significant unless one has a good 
quality spectrum with adequate signal-to-noise. The smaller 
one makes the spatial resolution element (in order to probe 
finer and finer regions), the larger the collecting area re- 
quired in order to obtain good-quality spectra in reasonable 
time-scales. Good spectral resolution will also be necessary 
to separate the effects of the varying enrichments and allow 
one to constrain two-temperature models with freely vary- 
ing abundances for each component. This will depend on just 
how extreme the metal variations are. Given that this is all 
for a regular, one-dimensional system, with no background 
file added, one can begin to appreciate the difficulties that 
genuine observations present in this regard. 

Accepting that direct detection of small-scale ICM 
metallicity deviations will be so difficult, we are forced to 
consider what forms of circumstantial evidence it may prove 
possible to employ. 

4-1.2 Abundance gradients 

Fig. illustrates that a bimodal abundance distribution 
within a cooling flow scenario has the ability to produce an 
apparent metallicity gradient where one does not in reality 
exist (the real mean abundance in our simulations remains 
unchanged, at least until metal-rich gas begins to drop out 
in the inner regions). For example, consider the results after 
1.0 Gyr of evolution, when the central abundance has risen 
to around 0.5 Zq. At later times the gradient effect would 
be even more pronounced if one were not able to resolve the 
central abundance drop. 

Large-scale abundance gradients were detected in sev- 



1997); A496 



( Matsumoto et al. 1996); AWM 7 (Ezawaetal 
(Dupke & White III 2000). The presence of an abundance 
gradient appears to be correlated with the existence of a 
cooling flow (e.g., De Grandi & Molend: [2001 ), although this 
may not be a true correlation (cooling flows generate gradi- 
ents), but merely due to the mergers which disrupt cooling 
flows also erasing any abundance gradient. Here we have 
demonstrated one mechanism by which cooling flows may 
actually give rise to the appearance of such abundance gra- 
dients. The interplay between cooling flows and abundance 
gradients has been studied previou sly by several authors. 
For example, Allen fc Fabian (|l99§l ) argued that it was the 
presence of abundance gradients in cooling flow systems that 
leads to the higher emission-weighted metallicitie s in these 



sources a s comp ared to non-cooling flow clusters. Reiseneg- 



ger et al. (1996) investigated the ability of cooling flows to 



create metallicity gradients by transport of metals. 

Recently, Chandra studies have confirmed the presence 
of large-scale abundance gradients in many clusters. The un- 
precedented spatial resolution of Chandra has revealed more 
detail in the profiles. In several cases, clusters are found 
to exhibit 'peaked' abundance profiles (i.e. a positive radial 
abundance gradient in the innermost regions, coupled with 
a negative gradient further out), very similar to th ose pro- 
duced in our simulations. For example: Centaurus ( [Sanders 



& Fabian 2002), where the negative abundance gradient 
peaks at 1.3-1.8 Zq at a radiu s ~ 15 kpc before falli ng back 



to 0.4 Z at the centre; A2199 (Johnstone et al. 2002), where 
the metallicity rises from ~ 0.3 Zq at 200 kpc to ~ 0.7 Zq 
at 30 kpc before dropping bac k to 0.3 Zq w ithin the central 



5 kpc; and possibly Perseus (Schmidt ct al. 2002), where 
there may be a high metallicity ring of around O.6Z0at a 
radius of 60 kpc. It is of course intriguing when a hypothesis 
designed to answer one issue (lack of low temperature line 
emission) ends up providing a potential explanation for an- 
other (abundance gradients with off-centre peaks); though 
of course there is no shortage of alternative explanations for 
both these phenomena. 

Extrapolating our results for several Gyr, one would 
predict that over not too long a timescale, the central region 
would become devoid of metals, with a ring of high metallic- 
ity at large (hence easily observable) radii that would have 
been detected before the Chandra era. This is clearly un- 
physical, and is easily explained away by our simple treat- 
ment of a complex problem. We start with fully formed 
metallicity variations rather than allowing them to develop; 
the variation is extreme - rather than allowing for several 
'phases' of differing abundances we use only two, one with 
no metals and one with a high abundance; and we do not 
account for replenishment of metals through continued en- 
richment, which would offset the extremes of behaviour that 
our simple models predict. 

The metallicity distribution in Perseus a ppears to show 



no cor relation with the galaxies. Intriguingly, Schmidt et al , 
( [2002[ ) also report what may be the first signs of small-scale 
metallicity variations in the Perseus ICM, although the scale 
of the effect is at the limit of detectability with current ex- 
posure times and is not statistically significant. 

A similar off-centre abundance peak is seen in the Chan 



dra observation of the merging cluster A3266 ( Henriksen 
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Figure 8. The temperature dependence of the equivalent widths 
for various ICM plasma spectral lines, for a solar metallicity 
plasma. Dot-dashed iron K; dashed iron L; dotted oxygen K; solid 
silicon K. See text for details. 



Figure 9. Ratios of element abundance to that of iron. Results 
are for the last of the times plotted in Fig. ^. Filled circles - 
magnesium; open squares - neon; open circles - oxygen; open 
diamonds - silicon. 



& Tittley 2005). These authors suggest an alternative ex- 
planation for this feature involving the merging subcluster 
depositing enriched gas close to the centre of the cluster. 
Shock heating is then called upon to preferentially devoid 
the enriched region of the lightest, most mobile elements, H 
and He, so that the metallicity may be increased without 
an associated raising of the density, for which there is no 
observational evidence. This explanation requires efficient 
motion of ions, something which it is not clear can happen 
even if conduction is relatively unimpeded. Moreover, these 
efficient transport properties are required precisely in the re- 
gions where subcluster merging is taking place. Cold-fronts 
(see Section ||) are plausible evidence that transport proper- 
ties may be highly suppressed in such volumes, presumably 
due to separate magnetic structures existing in the merg- 
ing components. It is intriguing that off-centre abundance 
peaks appear to be found both in clusters with and without 
cooling flows. If we are looking at the same phenomenon in 
both cases (this of course by no means clear), and if a com- 
mon mechanism is responsible, this would suggest that it is 
something associated with neither the cooling flow process, 
nor the merger process. 



4-1.3 Equivalent width effects 

The original A SCA abundance gradients were clearly de- 
tected in the equ ivalent width of the i ron K line (e.g. 
Centaurus where 



Fukazawa et al 



1994 report an increase 



by a factor of three on moving to the central regions of the 
cluster). Excluding the iron L complex from the Chandra 
data for Centaurus by simply fitting to the high-energy end 
of the spectrum (3-7 keV), whilst obviously increasing the 
noise, does not affect the qualitative nature of the conclu- 



sions regarding the abundance profile (Sanders & Fabian 
200^ ). Fitting our model spectra just using the data in this 



region produces different results, however, as illustrated in 
Fig. ^3. The noise level has increased somewhat, but the 
qualitative behaviour has also altered. There is still some- 
thing of an initial increase in the central abundance at early 
times, leading once again to the generation of an apparent 
gradient. At later times, however, the highly peaked abun- 
dance profile that was present when fitting to the entire 



spectrum does not appear. Instead the central abundance 
merely dies away. Excluding the iron K lines, however, by 
fitting only to the spectra between 0.3-5 keV produces much 
the same results as using the entire spectrum. 

Thus there is a qualitative difference in behaviour when 
fitting to the iron K lines as compared to the iron L lines that 
must be explained. In Fig. ^ we illustrate the temperature 
dependence of the equivalent width for various important 
ICM spectral lines, calculated from mekal spectra. 

Consider just the results for the iron. Upon cooling from 
high temperatures, the equivalent width of the iron K lines 
increases somewhat, then rapidly dies away below 1 keV. 
The width of the L complex, on the other hand, increases 
strongly and monotonically as the temperature decreases. 
Consequently, with a bimodal metallicity, the strength of 
the iron K lines relative to the continuum increases to a rel- 
atively small extent on cooling, producing the initial slight 
increase in the central abundance, before dying away as cool- 
ing progresses. The strong monotonic increase in the iron L 
width with cooling produces a much greater change relative 
to the continuum. 

To illustrate this point further, the equivalent widths for 
the silicon and oxygen K lines are also plotted in Fig. ^| Both 
increase monotonically as the temperature is reduced, but 
the rate of change for oxygen is relatively small, whereas sil- 
icon shows very similar behaviour to iron L. The arguments 
of the previous paragraph are borne out by the results of 
vmekal fits to the model spectra (Figs. |^c and d show the 
fits for silicon and oxygen respectively). The qualitative be- 
haviour of both profiles is the same, but the oxygen profile 
shows a less pronounced peak. 

In Fig. bl we plot the ratio of element abundance to iron 
abundance (both in solar units), for various elements, for the 
last of the times illustrated in Fig. ^. Neon and magnesium 
(and to a lesser extent silicon) appear overabundant, whilst 
oxygen appears slightly underabundant. 



4.2 Genuine abundance gradients 

Given that flat metallicity profiles coupled with bimodal dis- 
tributions are unable to fully reproduce the observed abun- 




Figure 10. Time evolution of the observed abundance profile for clusters with genuine abundance gradients. Vertical error bars are ltr, 
simulated observation time 25 ks, all results for single temperature mekal models unless otherwise stated. Top two panels, bimodal metal 
distribution: filled circles - 0.0 Gyr; open squares — 0.4; open circles — 0.5; open diamonds - 0.8. Bottom two panels, homogeneous metal 
distribution: filled circles — 0.0 Gyr; open squares — 1.3; open circles — 2.5; open diamonds — 3.8. Left-hand panels: fits to 0.3— 7.0keV 
region, dotted lines show results from two temperature fits where the reduced \ 2 f° r the single temperature fits exceeded 2.0. Right-hand 
panels: fits to the 3.0-7.0 keV (i.e. iron K) spectral region. 



dance trends, we have considered also clusters with genuine 
gradients in their abundance profiles. We have parameter- 
ized the profile in the following manner: 

z(r) = J Z+(Z c ~Z)(l~t-J r<r (n>0) 
\ Z r ^ r 

where Z = 0.3 is the mean cluster abundance; Z c = 2.0 
is the central abundance; ro = 200 kpc is the extent of the 
gradient; and n = 1.5 controls its severity. 

We have carried out simulations both with homoge- 
neous metallicities (i.e. the distribution smoothly follows the 
profile), and with a bimodal behaviour similar to before (i.e. 
90 per cent of the gas is metal-poor, the other 10 per cent has 
ten times the abundance dictated by Eqn. |l^). The results 
are displayed in Fig. |l^. Once again, the bimodal system de- 
velops an off-centre abundance peak, of a larger amplitude 
than before. This time, the trend is reproduced when fitting 
solely to the high-energy part of the range, panel (b) . In the 
homogeneous system, on the other hand, no such peak de- 
velops, rather the abundance declines reasonably smoothly 
back down towards a flatter state. Note that the temporal 



evolution is in both cases faster, owing to the higher abun- 
dances that are present. 

4.3 Cooling flow equivalent widths 

We have used simulations similar to those described in the 
previous section to produce synthetic XMM-Newton RGS 
spectra in order to investigate in more detail the ideas of 
Section ^| on the suppression of the low-temperature lines. 
These broadly reproduce the same results as the isobaric 
model, Eqn. ^| They do reveal another factor, however, 
which conspires to reduce the line suppression. When ge- 
ometrical (i.e. position within the cluster potential) consid- 
erations are taken into account, the metal-rich gas always 
has a larger cooling radius at any instant than the metal- 
poor gas, as a fairly obvious consequence of its faster cooling 
rate. In the case of an 8.6 keV cluster, for example, where 10 
per cent of the gas is made to be metal-rich, the fraction of 
cooled mass originating from the enriched phase asymptotes 
to about 1/8, rather than 1/10. This reduces still further the 
equivalent width suppression we can expect from a bimodal 
distribution of metals. 
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5 DISCUSSION 

To address the question of how small-scale ICM abundance 
variations might come about, supernovae can be considered 
as essentially point sources (on ICM scales) of extremely 
high metallicity gas. In of the mselves, these have structures 
that are far from simple (e.g., Hughes et al. 200C| ). One way 
in which escape of the enriched gas from the potential well 
of the host galaxy in to the ICM may occur is via superwinds 



(e.g., Heckman 2001) (although other processes such as ram- 
pressure stripping and pregalactic winds are doubtless signif- 
icant). Superwinds have a complex, multiphase structure in 
which it is difficult to probe directly the energetic, enriched 
gas driving the wind. There is, however, clear evidence for 
complicated structure in the wind on very sma ll scales in the 
optical and soft X-ray (e.g., [Strickland 2001 ). In our opin- 
ion, the question ought to be posed the other way around: 
how might the ICM become uniformly enriched? 

There are essentially two issues that determine whether 
or not the situation as we have chosen to model it is physi- 
cally realistic: i) how are the metals injected into the ICM, 
and ii) once in the ICM, do the metals move to spread out 
over a wide area or remain confined? The answers to both 
these questions remain uncertai n at pre s ent. T o address the 
second question first, following [Spitzei (1962), we may ex- 
press the deflection time (average time for a cumulative de- 
flection of 90°) for particles of mass m, charge Ze diffusing 
through field particles of mass mi, charge Zfe, density m as 



mate that iron nuclei may diffuse about 20 kpc in 1 Gyr. To 
maintain small-scale abundance variations over significant 
timescales in the ICM, we therefore require a strong sup- 
pression of diffusion. This result is to be expected, as it is 
equivalent to the strong suppression of conduction that has 
traditionally been invoked for the multiphase cooling flow 
picture. Magnetic fields have normally been appealed to as 
the causative agent to dampen tra nsport prope rties in an 
ionized gas. Early calculations (e.g., Tribblej 1989| ) appeared 
to show that a tangled magnetic field would indeed produce 
a strong reduction in transport. More recent calculations 
considering a chaotic field with turbulence extending over 
a wide range of length s cales indicate that the supp ression 
may in fact be minimal ( Narayan fe Medvedcv 2001 ). 
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with erf the standard error function, In A the Coulomb log- 
arithm, and assuming that all species are in thermal equi- 
librium at temperature T. The factor T ^ 1 is the Chan- 



In practice, simple diffusion is probably unlikely to be 
the limiting factor for the spread of ICM metals. Convec- 
tion, turbulence due to radio sources, galactic wakes, etc., 
will also play roles to varying degrees. If, however, the abun- 
dance drops observed with Chandra in the central regions 
of the ICM for several clusters are genuine, then this will 
imply limits on the amount of convection or mixing that 
can have taken place (else these features would have been 
smoothed out). It does not seem unreasonable to explore 
some of the consequences of small-scale metallicity varia- 
tions, and to keep the possibility in mind during spectral 
analyses. 

The question as to whether or not the metals in the 
intracluster and intergalactic media are homogeneously dis- 
tributed remains an open one. There is a non-negligible 
scatter in Galactic stellar metallicities for stars of all 
ages (the scatter increases for low metallicity stars), which 
has been taken to suggest that that Galactic disc has 
been chemically inhom ogeneous throughout its development 
(e.g., McWilliam 1997 ). Classical galactic chemical evolution 
models have tended to assume instantaneous dispersion of 
synthesized elements. 

The average metallicity of the intracluster medium on 
large scales has for some time been reasonably well estab- 
l ished at roughly 1/3 solar in most cases, both for nearby 
(|Edge fe Stewart 1991) and distant (Mushotzky & Loewen- 
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drasekt jar correction tor the hmte mass ot the held particles. stein j | 1997 p caters. Until comparatively recently, however, 
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with A the test particl e mass in units of the proton mass. 

As was shown by Rephaeli ( 1978] ) for the case of sed- 
imentation, in the ICM the contributions from field par- 
ticles other than protons (particularly helium nuclei) are 
significant. Summing over the elements in a 0.3 Zq plasma 
with iron nuclei as test particles, we find that r is reduced 
by around a factor of two from the value for a pure pro- 



there has been little information on how the ICM metals 
might be distributed on finer scales. There is, in our opin- 
ion, no particular reason why the whole ICM should be uni- 
formly enriched to the same metallicity, although this might 
be one's natural assumption. As we discuss in Section 4.1.1, 



direct detection of small-scale abundance variations in the 
ICM will be very difficult. Consequently, we feel it is impos- 
sible to rule out such variations at the present time. 

Recently, a deal of support has been given to the idea 
that thermal conduction might be operating in cluste r cores 



ton plafcma. At 10 K, and with nn — 10 cm , we hnd 
the corrected deflection times for iron and helium nuclei are 
T Fc » 0.3 Myr, r Ho « 4 Myr. 

From standard kinetic theory, the root-mean-square 
three-dimensional distance a particle will diffuse in a time t 
is given by 



at significant levels (e.g., |Narayan fe Medvedev 



et al 



2002| ; [Fabian et alj|2002| ). Whilst this runs c ontrary 



2001 



Voigt 



to the established multiphase cooling flow picture (Nulsen 



198f), it seems to have a degree of success in matching 



the observed temperature profiles. It is certainly the case 
though, that at least in some regions of the ICM, conduc- 
tion is highly suppressed. This is revealed by the 'cold-fronts 



6k B Trt 



(18) 



seen in several cluster c ores, e.g. A2142 ( Markevitch et al 
2000j). It was shown by |Ettori fe Fabian fl2000|) that such 



where A = Tv Ilns is the mean free path. Hence we can esti- 



features, which correspond to a sharp jump in the surface 
brightness profile, indicate a strong (factors of several hun- 
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dred) local suppression of the thermal conductivity, else the 
temperature discontinuity that is responsible for the bright- 
ness jump would be quickly washed out. If thermal conduc- 
tion is suppressed, then ion movement must be even more re- 
stricted, since net heat conduction may still take place with- 
out individual electrons travelling a significant distance. If 
thermal conduction is operating with a degree of efficiency, 
however, then the same need not necessarily apply to ion 
motion. If this were the case, though, there would be the 
possibility of sedimentation of the heavy elements (Fabian 



& Prin@ |l977| ; |Rephaeli| [l978 



Fabian et al 



|2002|) . Cold- 



fronts of course only provide direct evidence for inhibited 
thermal conduction across the fronts themselves. It is pos- 
sible (tho ugh as yet far from proven) that con duction may 



reach the spitzer (1962) value elsewhere (e.g., Voigt et al 



2002). If conduction is operating efficiently, it does not nec- 
essarily rule out small-scale metallicity variations, but it 
does weaken the case, both from the point of view of re- 
ducing/removing the motivation for them (suppressing the 
low-temperature cooling flow line emission by suppressing 
the low temperature gas), and by making it less likely that 
such conditions can persist. 

Our modelling is simplistic, and it may be argued that 
the conditions as we have chosen to represent them are not 
physically relevant. For example, there is no possibility of 
segregation of elements based on weight (that is, it is not 
possible for the metal- rich gas to sink to the centre). This is 
a consequenc e of our adopting the theoretical framework of 
Nulsen (1986), which requires that in a multiphase flow the 
various phases co-move, i.e. there is a single velocity profile. 



The a rguments for this idea are set out in detail in Nulsen 
(1986). Also, the timescale over which the changes in the 



abundance profile occur is somewhat short. This is a result 
of the rather unphysical initial condition; starting out with 
regions enriched to 3 Zq and not allowing any replacement 
of cooling metals. A more realistic treatment would allow 
the metallicity to build up initially with time and then al- 
low for some replenishment, and possibly also for different 
distributions of SNe la and II products. However, we are 
not concerned here with producing detailed models of the 
evolution, but rather in seeing what the general trends of 
behaviour might be. Note that because we have restricted 
our attention to single-phase models, there is no real trans- 
port of material by the cooling flow, since all the mass loss 
(except for the m e tal-ric h zones) takes place in the centre. 
Reisenegger et al. (1996) demonstrated the ability of multi- 
phase cooling flows to create (genuine) abundance gradients 
through transport of injected materials. Thus we may ex- 
pect that a more sophisticated multiphase treatment of this 
process might reveal more complex effects. 

Given that the abundance gradients observed with 
ASCA and Chan dra ar e seen clearly in the iron K line, the 
results of Section 4.1.; imply that this cannot be solely due 
to a bimodal distribution of metallicities. The results are 
not, however, inconsistent with such a metal distribution, at 
least in the simple scenario outlined here. At the very least, 
we have presented another reason to be wary of the iron 
L co mplex when fitting to X-ray spectra (e.g., Finoguenov 



et al. 200C). In our case, it is not due to the complex atomic 
physics of the L shell (the code we use to generate the spec- 
tra is the same as that we use to fit them), but rather to 
the temperature dependence of the equivalent width. Recall 



that iron is really the only ICM element for which one has 
two strong spectral indicators (K and L) ; for other elements 
only the K lines are useful. Playing devil's advocate, one 
could therefore imagine a situation in which the observed 
gradients in the iron profile are genuine, but those in the 
other elements are due to a process such as the one outlined 
here. This would, for example, severely impact on estimates 
of the SNe Ia:II importance ratio. Of course, this simple 
model would not explain the matching radial scales of the 
variations for iron and the other elements in this case. Nor 
could it explain any correlation between the iron profile and 
the visible light in the central regions of the cluster. 

Another process which may give rise to the same sort 
of radial abundance profile as those seen here (namely a 
peak in the abundance at an off- centre position) is resonant 
scattering (Gil'fanov et al. 1987). This acts to redistribute 
photons from the central regions of the cluster (where the 
o ptical depth is h i ghest) to a surrounding ring. See the work 
of 



Mathews et al 



001) for an application of these ideas to 
M87. Computation of detailed optical depths requires knowl- 
edge of the velocity structure of the gas along the line of 
sight. We will not comment on the resonant scattering is- 
sue here, except to say that to some extent the ideas of this 
paper and those of resonant scatte r ing a re in conflict. As 
was pointed out by Wise & Sarazin (1992), clumping of the 



ICM reduces the amount of scattering that takes place. X- 
ray surface brightness profiles depend on the rms density 
yj (n 2 ), whereas electron scattering depends on the mean 
density (n). The rms of a set of numbers is always greater 
than their mean (unless the numbers are all equal). Thus 
increasing the degree of clumping in the plasma reduces the 
mean density relative to the rms density and so reduces the 
relative effect of scattering. 

Another possibility is that the 'extra' metals in the 
central regions of man y clusters are to b e asso ciated with 
the central cD galaxy. Makishima et al. (2001) have sug- 



gested an alternative explanation for the enhanced emission 
seen at the centre of many clusters, which has traditionally 
been attributed to the cooling flow phenomenon by many 
researchers. Instead, these authors suggest the excess may 
be associated with the interstellar medium of the central 
cD galaxy. One argument invoked against the cooling flow 
interpretation is the fairly frequent presence of metallicity 
increases near cluster centres. Our present work suggests a 
mechanism by which such effects may indeed be produced 
by cooling flows. The reality or otherwise of any central dips 
in abundance would be an important discriminant for these 
two interpretations. 



Bohringer et al. ( 2002 ) also give some consideration to 



the possible effects of inhomogeneous metallicities on the ob- 
served abundances (section 2.2 of their paper). Unlike our 
consideration of individual emission lines, they look at the 
overall observed metallicity of the spectrum, as inferred from 
the global ratio of power in all emission lines to power in 
the continuum. These authors also discuss some compar- 
isons between the shape of the spectrum around the 1 keV 
region for M87 and the predictions of the bimodal model. 
All hypotheses live or die through comparison with data, so 
such investigations are highly necessary, but there is a large 
parameter space to investigate. And as they point out, such 
checks must be made on a source- by-source basis. They find 
a poor fit between the actual spectral shape and the pre- 
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dictions of the bimodal hypothesis (e.g. their fig. 6). Note, 
however, that this is for the case where the normalizations 
of the metal-rich and metal-poor phases are 'roughly equal'. 
We would not expect such a division to be successful in 
reducing the EW of the low temperature lines by an appre- 
ciable amount. 

Our main result is that the observed equivalent width 
suppression for low temperature emission lines in a cooling 
flow spectrum due to a metal distribution which is inhomo- 
geneous on small-scales is not as great as one might expect. 
For example, if all the metals are concentrated in 10% of the 
gas, th e suppression of the low-temperature lines relative to 



those fi om high temperature gas is only about a factor of 3 



(Section 2.2.2). It seems unlikely, therefore, that this method 



in isolation could produce a reduction in equivalent width 
equal to that seen in data, without pushing the bimodality 
of the metal distribution to extreme levels. There is an ef- 
fect, but it is not large enough. We have also shown that 
small-scale metallicity variations can give rise to in terest ing 
effects in the observed abundance profiles (Section 4.1.2) as 
compared to the true profiles. Such effects would give rise 
to serious difficulties in terms of interpreting abundances in 
cluster central regions. The possibility of small-scale metal- 
licity variations ought to be borne in mind when analysing 
high resolution X-ray spectra of cluster central regions. 
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